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Abstract 

This work deals with the possible solution of the solar neutrino problem in the 
framework of the resonant neutrino spin-flavor precession scenario. The event rate 
results from the solar neutrino experiments as well as the recoil electron energy 
spectrum from SuperKamiokande are used to constrain the free parameters of the 
neutrino in this model (Am^ and fi^) . We consider two kinds of magnetic profiles 
inside the sun. For both cases, a static and a twisting field are discussed. 



1 Introduction 

The amount of accurate solar neutrino data available at present, the numerous cheks of 
the functioning of the solar neutrino detectors that have been and are being performed 
together with more precise results in the field of solar modeling which are in a very 
impressive agreement with high-accuracy helioseismological data, suggest strongly that 
the observed deficiency of the solar neutrino is one of the most convincing indications 
for new physics beyond the standard electroweak theory. 

Two attractive solutions to this puzzle are (i) the Mikheyev - Smirnov - Wolfen- 
stein (MSW) [1] matter-enhanced neutrino oscillations, and (ii) the spin precession via 
a large magnetic moment of the electron neutrino(/i,y) proposed by Okun, Voloshin and 
Vysotsky (0VV)[2] motivated by the apparent anticorrelation of neutrino flux in the 
Davis experiment with the sunspot activity [3]. This precession is enhanced resonantly in 
the presence of matter yelding to a resonant spin flavor precession (RSFP)[4] analogous 
to the MSW effect. 
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The present paper deals with the RSFP type solution to the solar neutrino problem. 
This solution has been advocated by many authors[5,6]. The magnitude and the profile 
of the sun magnetic field beeing the most uncertain variables, they have been assigned 
a variety of values and forms (for a review see [5, 6]). It is worth noticing that the solar 
neutrino experimental results may be used to probe the solar interior (its magnetic field 
profile for instance) [7]. We explore the preliminary 708 days^ of operation results from 
the SuperKamiokande [8,9] experiment together with the last results of the four solar 
neutrino experiments (Homestake[10], Kamiokande[ll], SAGE[12] and GALLEX[13]) for 
two variants of the solar magnetic field. The range for the RSFP neutrino parameters 
Am^ and /ii, compatible with the total measured event rates in all of the solar neutrino 
experiments, as well as with the SuperKamiokande electron recoil energy spectrum are 
investigated both for a static and a twisting magnetic field. The paper is organised as 
follows: In section 2 a brief review of the experimental situation is presented. In section 
3 the neutrino propagation equation through solar matter is solved analytically using 
the Landau Zener[16] approximation. Hence we get the neutrino parameters allowed by 
the total event rates in neutrino experiments and the distortion of the recoil electron en- 
ergy spectrum measured by SuperKamiokande. The evolution equation and the allowed 
regions for neutrino parameters are next considered in section 4 taking into account the 
twisting effect of the sun magnetic field. Our conclusions are summarised in section 5. 
Some details about the contour finding procedure are presented in the appendix. 



2 Experimental status 

Let us first briefiy recall the experimental situation. Table 1 summarize the neutrino 
event rates that have been measured in the four pioneering solar neutrino experiments 
Homestake[10] , Kamiokande[ll], SAGE[12] and GALLEX[13] together with the neu- 
trino flux measured after 708 runing days by SuperKamiokande[8,9]. The observed event 
rates are significantly smaller than the theoretical expectation of the BP98 model[17]. 
The GALLEX and SAGE experiments measure the same quantity, in what follows we 
consider only their weighted average rate. We also adopt the SuperKamiokande mea- 
surement as the most precise direct determination of the higher energy neutrino 
fiux. 



3 Static field 



3.1 Time evolution of the solar neutrino 



Disregarding, for simplicity, possible neutrino mixing, the time evolution equation of 
the neutrino in the transverse magnetic field B is expressed in the Majorana^ weak 
interaction (I'ci^n) sector as: 
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^The present statistics of SuperKamiokande is about 825 days. 
^For Dirac neutrino, just replace | in the Gp-ne coefficient by 



where Am^ is the neutrino flavour mass square difference, and Gt- is the Fermi couphng 
constant. We used the approximation n„ ~ ^rig between the electron and neutron den- 
sities in the sun, vahd in the convection and upper radiation zones where most of the 
neutrino trajectory hes[14]. 

The electron density Ue decreases exponentially along the neutrino trajectory and is well 
approximated by[14]. 

GF.ne = 2.11xlO-"exp(-^^^) eV, (2) 

r being the distance from the center of the sun and Rq being the solar radius. 
We use the analysis of Parke[15] to determine the average survival probability: 

P^e-^u, ^l + i^- Plz) cos 2§i cos 20 f, (3) 
where denotes the initial(final) value of the mixing angle 9 such that, 

sin^ 9 = '-^ 2 > (4) 

and the jump probability P^z is taken in the Landau Zener'^ approximation[16] 
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calculated at the resonance point obtained by requiring, 
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The generalisation of eq.3 to include the neutrino energy and production range dis- 
tributions is thus: 

f 1 1 . . 

P^,-.u, = J WEWri- + ( 2 - Plz)) cos 29i cos 29fdEdri, (7) 

where the function we, Wr represent respectively the probability density of neutrino 
production per unit energy and per unit length. 

Unfortunately, the magnetic fleld inside the sun is not accessible to direct observation. 
At the moment there is no model for the solar magnetic fleld and very little is known 
about it: not only its proflle is unknown, but even its strenght is very uncertain. 
An upper limit on the strengh of the solar magnetic field comes from the requirement 
that the field pressure must be smaller than that of matter. This is a rather weak limit 
(for the convective zone B < lO^G) and all other more stringent bounds arc highly model 
dependent. Several plausible profiles have been proposed and investigated in literature. 
In the following we consider two distributions of the solar magnetic field: 



^The Landau-Zener approximation which assumes a linearly decreasing density in the vicinity of the 
critical point works rather well in the sun[18] 



1. Linear distribution: 

= lO^G r < 0.7Rq , . 

= 105^(1 - 3.33(r/i?0 - 0.7)) r > 0.7i?o- 



2. Wood-Saxon (WS) distribution: 



B(r) = (9) 

^ ^ l + explO(r-i?o)/i?o' 



3.2 Total event rates 

Using only the total event rates measured at the Homestake, Gallium and Super- 
Kamiokande^ experiments, figure 1(2) shows the allowed region in the (/i;y,Am^) pa- 
rameters space for the linear (Wood-Saxon) profile. The black dot within each allowed 
region indicates the position of the best fit point in the parameters space. The best fit 
for the linear profile is obtained for: 



Am^ = 1.8 X IQ-^eV'^, 
Id = 3.9 X IO-^Vb, 



(10) 



for which Xmin — l-S*^- Wood-Saxon profile, the best fit occurs at: 

Am^ = 1.7 X 10-^2, 



II = 7.7 X IO-^Vb, ^^'^^ 



with xLn = 0.94. 



3.3 Recoil electron energy spectrum 

Unlike the neutrino event rate deficit, the energy spectrum of recoil electron observed 
at SuperKamiokande is one of the most important model independent solar observables. 
Therefore, a deviation of the observed electron recoil energy spectrum shape from what 
is predicted by standard electroweak theory would be an indication of new physics (such 
as neutrino oscillations). 

The SuperKamiokande experiment [9] has measured the energy spcctrTim of recoil 
electrons from the neutrino-electron elastic scattering in water above 5.5 MeV. In what 
follows we will use the SuperKamiokande data relative to 708 days of operation. The 
data are given in an 18-bin energy histogram, where for each bin we have the ratio 
between the experimental event rate and the theoretical one. The first 17 bins have a 
width of 0.5 MeV starting from 5.5 MeV while the last bin include events with energies 
from 14 MeV to 20 MeV. 

Fig. 3 and fig. 4 show (respectively for linear and W-S magnetic field profiles) the RSFP 
neutrino parameters regions that are allowed when we take into account the information 
from the SuperKamiokande recoil electron energy spectrum alone. For both profiles, a 

^Since the quoted uncertainty in the Kamiokande rate is much larger than the uncertainty in the 
SuperKamiokande, the results are essentially unchanged if the rate from kamiokande is also considered. 



large region of RSFP parameters space is consistent with the data. The best fit is 
obtained with the hnear magnetic field at: 



Am^ =3.0x10-^2, 

PL = 3.1 X IO-^Vb, ^ ' 

with Xmin = 0-90. 

For the Wood-Saxon profile the best fit occurs at: 

Am2 = 1.0 X 10-^e\/^ , . 

11 = 3.7 X IO-^Vb, ^^"^^ 

and xlain = 1-08. 

Fig. 5 and fig. 6 show the regions compatible with both the event rates and the infor- 
mation from SuperKamiokande spectrum data. Even though the best fit solution consid- 
ering only the spectrum information (the dark point in figures 3 and 4) do not lie within 
the allowed regions by the analysis taking into account the event rate information, a 
large region in the neutrino parameters space (Am^, ^) are consistent with constraints 
from the total measured event rates as well as with those from the SuperKamiokande 
recoil electron energy spectrum. In fig. 7 we plot the recoil electron energy distribution 
divided by the standard prediction expected to be observed in the SuperKamiokande 
detector, using the best fit parameters found (eqs.lO and 11 respectively for the hnear 
and the W-S profiles) together with the data from SuperKamiokande[9]. 



4 Twisting magnetic field 



Several works proposed that the transverse component of the solar magnetic field may 
change its direction along the neutrino trajectory[19]. This can lead to new interest- 
ing phenomena in neutrino physics[20]. In this section we report on the observed deficit 
interpreted in terms of RSFP for a twisting magnetic field. We also discuss how such rota- 
tion would affect the recoil electron energy spectrum observed by the SuperKamiokande 
experiment. 

The equation for the flavor neutrino wave function describing the propagation of 
neutrino in matter with twisting transverse magnetic field can be written -for the case 
of interest- as: 



d 
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The angle 0(t) defines the direction of the magnetic field Bxit) in the orthogonal plane 
to the neutrino momentum and Bt = \ B^ {t) \ . 

The evolution equation looks like the one obtained in the static case (no twisting) 
with the addition of a quantity proportional to to the effective matter density . The 
formalism sketched above for the calculus of the neutrino survival probability remains 
valid. 

Let us assume that (f) ^ ^ where rg is the curvature radius of the magnetic field lines, 
then it is easy to see that the effect of the twisting field becomes significant when: 

1 ^Gf 



For matter density values taken at the bottom of the convective zone of the sun, this 
gives: 

ro ~ O.li?0. (16) 

We did the same calculation as before assuming both linear and W-S profiles for 
Bxir-) and taking for k = Rq/tq (a dimensionless factor which characterise the twisting 
velocity) the values ±10 according to the interesting feature expected from eq.l6. 
Fig. 8 shows the allowed region obtained using the hnear profile while fig. 9 gives the 
analogous result for the W-S profile. 

For k=+10, the two field configurations used in our calculations give poor fits to the 
total event rates. (Xmm ~ ^-^"^ Xmin ~ ^-^^ the W-S and the linear profiles 
respectively) 

In contrast, when the field twists in the opposite side (k=-10), good fits to the total 
event rates are obtained with both configurations of the magnetic field. 
Using the hnear field distribution, the best fit is found for: 

Am^ = 9.9 X lO-i^ey^ 

pi = 2.3 X lO-^Vs, ^ ^ 

with a shallow Xmin — 0.18. In an analogous way, we find the best fit to the data for the 
W-S field at: 

Am^ = 8.4 X 10-l^ey^ , . 

li = 3.4 X IO-^Vb, ^ ' 

with xlin = 0.2. 

The most important change in the allowed range of neutrino parameters compared 
to the standard case (no twisting) is the disappearence of the area at the lower right 
corner of figs. 1-2 compatible with the measured experimental event rates at the 99%CL. 
The allowd regions are also extended to somewhat smaller values of Am^. The main 
reason for this is the fact that the resonant density depends on the magnitude of Am^ 
as well as on the velocity 4> which has a tendency to move it inward the sun, so for 
Am^ small enough, different solar neutrinos types can also undergo resonant transition 
in both convective and radiative zones. 

Fig. 10 and fig. 11 show the regions of neutrino parameters compatible with the Su- 
perKamiokande spectrum of the recoil electron energy for the field distributions given 
by (8) and (9) and the two twisting magnitudes (k=±10). 

In the case k=10, better fits to the SuperKamiokande electron spectrum are obtained 
compared to the total event rates case. For a linear distribution we have: 

Am^ = 1.0 X lO-W^, 

H = 2.9 X IO-^Vb, (19) 
Xlin =0.68, 

while for a W-S distribution, the result is: 

Aw? = 8.0 X IQ-^eV^, 

PL = 3.8 X 10-^Vb, (20) 

_ n 01 
A.min u.OJ., 



Good fits are also obtained for tlie case (k=:-10) and tliey give: 



7.6 X IQ-^eV^, 

2.7 X (21) 
1.45, 

3.1 X lO'^eV^, 

4.3 X lO-^Vs, (22) 
1.52, 

for the W-S field. 

This results can be easily understood if we consider the dependence of neutrino sur- 
vival probability on the twisting magnitude 4> (see fig. 12). For a given Am^, at energies 
of the solar neutrino spectrum low enough, (j) is very small with respect to the neutrino 
oscillation coeficient It follows that the twisting effect is absent at this energy scale. 
The twisting effect becomes significant for the highest energy part of the spectrum and 
depends on the sign of 0. RSFP is amplified when (k=-10) and suppressed for the op- 
posite sign. This dependence is in favour of a positive twisting (k=10) to fit the electron 
spectrum measured by the SuperKamiokande experiment. 

Comparing the regions allowed by the total event rates shown in figs. 8 and 9 with 
those allowed by the energy spectrum shown in Figs. 10 and 11, it follows (Fig. 13-14) 
that for k=-|-10, none of the field distribution is able to explain the whole data from 
underground experiments. For the case k=-10, there is a large region in (Am^, ^) 
compatible with both constraints for the two fields. For this later case, and using the 
best fit parameters found (eqs. 21-22), we plot in figure 15 the expected RSFP distortion 
of the SuperKamiokande recoil energy spectrum. 

We stress that we have also studied two extreme cases corresponding to slow {k < 1) 
and fast(A; > 100) twisting. While no apreciable change from the standard case (no 
twisting) has been found for the case of slow twisting, the highest values of k decouples 
the (z/g, z7^) system and consequently the fast twisting scenario is surely unable to explain 
the deficit found by the underground experiments. 

5 Conclusion: 

The results from the SuperKamiokande experiment opens a new area in the solar neu- 
trino studies. Using the Landau-Zener formahsm, we have investigated the RSFP way 
for a solution to the solar puzzle in the light of the latest experimental results as well 
as the theoretical predictions. 

We have identified the allowed regions of neutrino parameters for either static and twist- 
ing magnetic field in the sun. This was done first by considering the constraints from 
total event rate results. The obtained subset of the parameters that are consistent with 
the total rates was then confronted to the SuperKamiokande electron recoil energy spec- 
trum to extract possible neutrino parameter regions compatible with the whole data set. 
We found that the RSFP scenario can account for both the observed deficiency of solar 
neutrino fiux and the measured SuperKamiokande spectrum. This is the case of the 



for the hnear field and: 
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Am^ 
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static field and of the rotating one with negative k. Although it provides a good fit to 
the spectrum data, solar twisting field with positive k seems to be disfavoured to explain 
the whole data available from underground experiments. 

We should stress that the quahty of the fits depends also on the chosen field configu- 
ration. Differences between the two cases as large as one order of magnitude in the /i 
values are found. 

We emphasize that the somewhat vagueness of the conclusions drawn by this work are 
related to the poor knowledge we have of the field distribution and -in the case of a 
rotating field- of the magnitude and the sign of the rotation velocity 0. 
It is crucial to have a good enough knowledge of the sun magnetic field parameters in 
order to draw an ultimate conclusion to the RSFP contribution to the solution of the 
solar neutrino puzzle. 

It is also of great importance to compare the RSFP with other oscillation scenarios (i.e. 
the MSW, the just so and the vacuum solutions) to set the more likely solution. The 
present energy spectrum data are uncertain enough (especialy in hight energy part of the 
spectrum) to inhibit any clear statement. However, the growing experimental accuracy 
will - in a near future- allow to validate/exclude more easily the various proposed mod- 
els. The possible time variation (anticorrelation with the Wolf cycle of the sun activity) 
that may be observed in the SK signal should also help making the comparison. 



Appendix 

Given the experimental rates Ri and their uncertainties cTj, the is defined as: 

Xrates = i^i,th " Ri,exp] [crij,tot\~^ [Rj,th - Rj,exp] , (23) 

where i,j label the experiment type (Ga, CI and SK). The expected rates and the error 
matrix are derivated following the method of [21] and 

^ij,tot — ^ij,exp + ^ij,thJ (24) 

^ij,exp ~ ^ij'^i,expO'j,exp, (25) 

denoting by (Xi^exp, the experimental error given in table 1. 

(Tij^th is the sum of two contributions: the one coming from uncertainties on cross sections 
(GS) and the one from the uncertainties on astrophysical parameters (AP) 

(^ij,th = (^ij,cs + ^ij,AP- (26) 

If we write a particular rate event as:i?j = J2 ^ij^Pj 
then we have 

8 QD QD 8 

k,l=l fc=l 

where Ri^. = Cik(/)k are the partial rates and the sum runs over the eight relevant neutrino 
fiuxes {i.e.pp,pep, hep,^ Be,^ B,^^ N,^^ O and ^^F). 



Similarly we obtain: 

8 

'^ij,AP ^ij H RikRjAln(f)kAln(f)i. (28) 

k,l=l 

The AlfKpi are calculated using the Bahcall's code exportrates.f [22] with some minor 
modifications. 

The correlation matrix used in this work is: 



Experiment 




Correlation matrix 


Ga. 


1.000 




CI. 


0.671 


1.000 


Super Kamiokande 


0.687 


0.964 1.000 



On the other hand, for the study of the observed SuperKamiokande energy spectrum, 
we use the following x^' 

X%ec = E [(^S,,tH - S,,e.A [WinPSj,tH - Sj,^.,] , (29) 

where Si^th is the predicted event rate for the i-th energy bin, Si^exp is the measured rate 
and /3 is a free parameter which normalizes the predicted solar neutrino flux to the 
measured flux. 

The input neutrino fluxes have been taken from [17]. Asymmetric errors have been 
conservatively symmetrized to the largest one. We use the improved neutrino cross 
section for each detector given in [23] and the neutrino spectra given in [24]. Allowed 
regions in the (^,Am^) plan are obtained by finding the minimum x'^ plotting 
contours of constant = Xmin'^^x'^ where A^^ = 4.61 for 90% C.L., 5.99 for 95% C.L. 
and 9.21 for 99% C.L. All results are given as function of the reduced x^.(i.e.x^=x^/ndof. 

) 

Majorana neutrinos with {i/g — > i/^) have been used through this work. The results for 
Dirac neutrinos being practically identical. 
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Experiment Ref. 



Data 




Homestake [10] 

Kamiokande [11] 

SAGE [12] 

GALLEX [13] 



2.56 ±0.16 ±0.16 7.7: 
2.80±0.19±0.33 5.1£ 



2.44 ± 0.051^7 5. 



77.5 ± 6.2lt? 



SuperKamiokande [8,9] 



Table 1: Neutrino event rates measured by solar neutrino experiments, and correspond- 
ing predictions from the BP98 solar model[17]. The quoted errors are at la. 




Figure 1: For the linear field we show the regions of parameters Am^ and obtained 
by a fit to the total event rates only both at 95%C.L. (doted erea) and 99%C.L.(solide 
line). The best fit is indicated by dark filled circle. 
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Figure 2: Same as figure 1 using the Wood Saxon profile for tlie solar magnetic field. 




Figure 3: Considering only the SuperKamiokande energy spectrum, allowed region for 
neutrino parameters at 99%C.L.(solide line) and 95%C.L. (dashed line) are shown for 
the linear profile. The best fit is indicated bu a filled circle. 
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Figure 5: For the linear field(eq.8) we compare the allowed neutrino parameters regions 
obtained separately by the total event rates (solide line) and the SuperKamiokandc 
energy spectrum (dashed line). The filled circle is the best fit to the SuperKamiokandc 
spectrum (eq.l2) while the open circle is the best fit to the total event rates(eq.lO). The 
comparaison is done at 99%C.L.. 
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Figure 6: Same as in figure 5 for tlie Wood Saxon filed. 




Figure 7: We plot the recoil energy spectra expected from RSFP scenario using the best 
fit parameters(eqs. 10-11) for the linear field (solide line) and Wood Saxon filed (doted 
line) divided by the SSM expectation. SuperKamiokande data are also shown with error 
bars representing the statistical and systematical errors added in quadrature. The data 
were directly read from M.B. Smy in [9]. 




Figure 8: Regions of the parameters Am^ and obtained by a fit to the total event 
rates only, using the field distribution given by eq.8 and the two values of k(see the 
text for definition). The best fit is indicated by dark circles, results are at 99%C.L. and 
95%C.L. (doted erea). 




Figure 9: Same as in figure 8 for the solar field given by eq.9 




Figure 10: Considering only the SuperKamiokande data spectrum, the allowed regions 
for neutrino parameters at 99%C.L. (solide line) and 95%C.L. (dashed line) are shown 
for the linear profile and two values of k. 




Figure 11: Same as in figure 10 using the Wood Saxon field. 




Figure 12: Neutrino survival probability computed using the linear(a) and the Wood 
Saxon(b). The neutrino parameters are those of the static best fit found in eqs. 10-11. 
The solid lines correspond to the static case. Doted and dash-doted hues refer to the 
rotated case respectively with k = -10 and k = -|-10. 




Figure 13: For the linear field(eq.8) and the two values of k, we compare the allowed 
neutrino parameters regions obtained by considering only constraints coming from the 
total event rates(figure 8) with those coming from the SuperKamiokande data spectrum 
(figure 10). 




Figure 14: For the Wood Saxon field (eq.9) and the two values of k, we compare the 
allowed neutrino parameters regions obtained by considering only constraints coming 
from the total event rates (figure 9) with those coming from the SuperKamiokande data 
spectrum (figure 11). 




Figure 15: As in figure 7, for the linear field (solide line) and Wood Saxon field (doted 
line) and for k=-10, we plot the recoil electron energy spectra from RSFP scenario using 
the best fit parameters (eq. 21 and eq.22) divided by the SSM prediction. 



